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Motivation 


•  Feedback  cycle  between  liquid  rocket  engine  (LRE)  combustion  chamber  pressure 
perturbations  and  unsteady  combustion12 

•  Large  amplitude  fluctuations  in  pressure  and  combustion  heat  release  rates  =^> 
combustion  instability 


Flow 

Heat  Release 

Perturbation 

Ei 

Oscillation 

Pressure 

Perturbation 


'Harrje,  D.T.,  and  Reardon,  F.H..  Scientific  and  Technical  Information  Office,  National  Aeronautics  and  Space 
Administration,  NASASP-194,  1972. 

2Schadow,  K.C.,  Gutmark,  E.,  Parr,  T.P.,  Parr,  D.M.,  Wilson,  K.J.,  and  Crump,  J.H..  19th  AIAA  Fluid 
Dynamics,  Plasma  Dynamics  and  Lasers  Conference,  AIAA  1987-1326 
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Objective 


•  Impose  external  acoustic  perturbations,  and  examine  the  response  and  stability 
characteristic  of  shear-coaxial  injector  flow  to  pressure  perturbation 


Acoustic/Pressure 

Perturbation 


Shear-Coaxial 
Injector  Flow 


•  Investigate  influence  of  injector  geometry  on  flow  response  to  external  pressure 
perturbation 

•  Vary  the  outer- to-inner  jet  momentum  flux  ratio,  /,  under  subcritical  and  nearcritical 
chamber  pressure  conditions,  i.e.,  reduced  pressures  Pr  =  0.44,  1.05 


Pr  = 


chamber 


critical,  N- 


Pcnuca,  n2  =493psi(3.4MPa) 


•  Apply  proper  orthogonal  decomposition  of  high-speed  image  pixel  intensity 

fluctuations  to  extract  spatial  and  temporal  characteristics  of  prevalent  coherent  flow 
structures 
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N  on-Dimensional 
Spatial  Growth  Rate 


Previous  Works  on  Jet  Instability 


Michalke  and  Hermann  (1982)  did  linear,  inviscid  instability  analysis  of  a  circular  jet 
with  coflow 

Showed  that  with  increasing  coflow  velocity,  Urx 

—  Helical  disturbances  more  unstable  than  axisymmetric  ones  farther  downstream  of  exit 

—  Jet  flow  becomes  less  unstable,  but  spectrum  of  spatial  growth  rate  becomes  broader  and  the  peak  shifts 
to  higher  frequencies 

Increasing  Downstream  Distance 


A 


Helical 


Axisymmetric 


->  Non-Dimensional  Frequency 


•  Dahm  et  al  (1992),  Wicker  and  Eaton  (1994)  conducted  experimental  investigation  of 
large-scale  vortex  structures  in  the  near  field  of  coaxial  jets 

—  For  outer-to-inner  jet  velocity  ratios  greater  than  one,  found  that  coherent  structures  in  the  outer  shear 
layer  dominate  those  in  the  inner  shear  layer 
—  At  large  axial  distances,  shear-layer  vortices  exhibit  helical  structures 


Dahm  et  al JFM  1992 
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Schematic  of  Experimental  Facility 
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Image  of  Experimental  Facility 


Piezo-Siren 


Thermocouple  Probe 


Waveguide 


Inner 

Chamber 


Differential  Pressure 
Transducers 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


High-Speed 

Camera 


Coaxial  Injector 


Injector  Configuration 


•  Two  types  of  outer- to-inner  jet  cross-sectional  area  ratios 
—  Large  Area  Ratio  (LAR) 

—  Small  Area  Ratio  (SAR) 


Injector 

t 

t/D1 

D} 
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2.44 

3.94 
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Acoustic  Field  Set-Up:  Pressure  Antinode 


•  Pressure  antinode  (PAN)  -  condition  of  maximum  pressure  perturbation  in  the 
acoustic  field 

•  Piezo-sirens  forced  in-phase 

•  Superposition  of  quasi- ID  acoustic  waves  traveling  in  opposite  directions  =>  PAN  at 
the  jet  location  (geometric  center  of  test  section) 


Proper  Orthogonal  Decomposition 


•  Proper  Orthogonal  Decomposition  (POD)  or  Principal  Component  Analysis 
(PCA)  was  used  for  extracting  dominant  dynamical  processes  embedded  in  high¬ 
speed  images. 

•  A  time-resolved  set  of  images  A(x,t)  can  be  represented  as  a  linear  combination 
of  orthonormal  basis  functions  (f)k  (aka  proper  orthogonal  modes)12  : 

M 

A(x,t)  =  ^ak(t)<f>k(x) 

k= 1 

where  ak(t)  are  time  dependent  orthonormal  amplitude  coefficients  and  M  is  the 
number  of  modes 


•  Main  idea:  POD  modal  amplitudes  capture  the  maximum  possible  “energy”  in 
an  average  sense3,  i.e., 

k  k 

where  bk(t)  are  the  temporal  coefficients  of  a  decomposition  with  respect  to  an 
arbitrary  orthonormal  basis  y/t 


1  Chatterjee,  A.  Current  Science ,  Vol.  78,  No.  7  (2000) 

2  Arienti,  M,  and  Soteriou,  M.C..  Phys.  Fluids  21,  112104  (2009) 

3  Narayanan, V.,  Lightfoot,  M.D.A,  Schumaker,  S.A.,  Danczyk,  S.A.,  and  Eilers,  B..  ILASS  Americas,  2011 

4  Berkooz,  G.,  Holmes,  P.,  and  Lumley,  J.L..  Annu.  Rev.  Fluid Mech.  25.  539  (1993) 
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Construction  of  Data  Set 


•  First,  form  a  row  vector  consisting  of  all  pixel  intensity  values  of  each 
snapshot  image  (with  resolution  of  n  rows  by  m  columns)  in  order  of  increasing 
columns,  then  increasing  rows 


Image  Frame 


columns 


I 

n 

rows 

i 


•  Then,  combine  all  such  row 
vectors  for  N  sequences  of  image 
frames  resulting  in  a  matrix  A 
consisting  of  N  rows  by  P  =  n  x  m 
columns  of  intensity  values. 
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Orthogonal  Decomposition  Technique 


•  Eigenvalue  decomposition  or  singular  value  decomposition  (SVD)  can  be  used 

•  SVD  preferred  since 

1.  Applicable  to  non-square  matrices  (most  likely  the  case) 

2.  Decomposition  matrices  are  orthogonal 

3.  Subroutine  readily  available  in  MATLAB® 

•  Subtracted  temporal  mean  of  A  =>  matrix  of  intensity  fluctuations  A 


•  Applied  SVD 


=>  A  =  US  VT  =  QVT 

Columns  of  Q  ~  ak(t)  contain  temporal  information 
Columns  of  V  ~  (f)k(x)  contain  spatial  information 


N  x  N  < 

Orthogonal  Matrix  of  Left 
Singular  Column  Vectors  of  A 


N  x  p 

Orthogonal  Matrix  of  Right  Singular 
Column  Vectors  of  A  proper 
orthogonal  modes  (POM) 


N  *  N 

Diagonal  Matrix  of 
Singular  Values 
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Results  -  Subcritical  Baseline  at  Low  J 


LAR,  Pr  =  0.44,  J  =  0.5 

Snapshot  Average 
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Power  Spectral 
Densities  (PSD) 
of  Temporal 
Coefficients  of 
POMs  1  and  2 


Results  -  Subcritical  PAN  at  Low  J 


•  LAR,  Pr  =  0.44,  /  =  0.5,  forcing  Frequency, /F  =  3.14  kHz 


Symmetric  Structures 

Identified  with  Characteristic 
Frequency  at  fF 


Cross-Power  Spectral  Density  (CPSD) 


•  CPSD  yields  the  FFT  of  the  cross-correlation  of  the  temporal  coefficients 

•  Magnitude  and  phase  plots  used  to  determine  existence  of  propagating  structures 


LAR,  Pr  =  0.44,  J  =  0.5 


PAN 
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Sample  Animation  -  PAN  (fF  =  3.14  kHz) 


y|jfcs  I  6MC0VW1  ■  C>tyUC#  i  Gfly Vtm 


•  LAR  Pr  =  0.44,  J  =  0.5 


Raw  Images  POMs  1  and  2  POM  1  POM  2 

Time  =  0.00  ms  Time  =  0.00  ms  Time  =  0.00  ms  Time  =  0.00  ms 


Superposition  of  POMs  1  and  2  Resulted  in  Downstream  Propagating  Structures 
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Results  -  LAR,  Pr  =  0.44,  Baseline 


•  Antisymmetric  flow  structures  indicated  helical  type  flow  instabilities  for  all  J 

j  =  2.1  J  =  5.2  J  =  11  J  =  20 


Results  -  LAR,  Pr  =  0.44,  PAN 


•  Gradual  shift  from  symmetric  to  antisymmetric  flow  structures  with  increasing  J 

•  Response  at  forcing  frequency, /F,  dominant  at  lower  J 


j  -  2.1  J  =  5.2  /  =  11  J  =  20 


Results  -  LAR,  Pr  =  1.05,  Baseline 


Antisymmetric  flow  structures  indicated  helical  type  flow  instabilities  for  all  7 
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Results  -  LAR,  Pr  =  1.05,  PAN 


•  Trend  in  response  with  varying  J  similar  to  Pr  =  0.44 

•  Gradual  shift  from  symmetric  to  antisymmetric  flow  structures  with  increasing  J 

•  Response  at/F  still  took  over  natural  (baseline)  frequency  at  lower  J 


Results  -  SAR,  Pr  =  0.44,  Baseline 


•  Helical  type  flow  instabilities  became  more  well-defined  with  increasing  J 
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Unlike  LAR  flows,  characteristic  peaks  showed  minimal  variation  in 
_ frequency  with  outer  jet  velocity _ 
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Results  -  SAR,  Pr  =  0.44,  PAN 


•  Symmetric  structures  persist  despite  increasing  7 

•  Response  at/F  strong  at  highest  7 


7  =  2.0  7  =  5.2  7=12  7=17 


Results  -  SAR,  Pr  =  1.05,  Baseline 


•  Antisymmetric  flow  structures  indicated  helical  type  flow  instabilities 


0Dr 


4  Dr 


8Z)jJ 


7  =  9.4  7=19 


Similar  to  Pr  =  0.44,  characteristic  peaks  showed  minimal  variation  in 
_ frequency  with  increasing  outer  jet  velocity _ 
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Results  -  SAR,  Pr  =  1.05,  PAN 


Similar  to  Pr  =  0.44,  symmetric  structures  persist  even  at  high  J 
Vortex-pairing  interactions  were  most  dominant  response  at  0.5/F 
Response  at/F  strong  at  highest  J 


j  =  9.4 

(fF  =  3.07  kHz) 


Conclusion 


•  Proper  orthogonal  decomposition  of  high-speed  image  intensity  fluctuation  data 
revealed  key  spatial  and  temporal  characteristics  of  flow  structures 

•  In  both  pressure  regimes,  LAR  injector: 

—  Peak  frequencies  of  baseline  flow  instabilities  became  broader  and  shifted  to 
higher  frequencies  with  increasing  J 

—  PAN  forcing  at  low  J  produced  symmetric  flow  structures,  while  at  higher  /, 
influence  of  forcing  subsided 

—  Spectral  magnitude  plots  showed  decreasing  influence  of  PAN  forcing  with 

increasing  J 

•  In  both  pressure  regimes,  SAR  injector: 

—  Increasing  J  had  minimal  influence  on  peak  frequencies  of  baseline  flow 
instabilities 

—  PAN  forcing  produced  symmetric  flow  structures  regardless  of  J 
—  Spectral  plots  showed  strong  response  to  PAN  forcing  at  low  and  high  J 

•  Operated  at  high  enough  /,  LAR  injector  flows  less  vulnerable  to  external  pressure 
disturbances 
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Data  Summary  Tables  -  LAR 


Pr  =  0.44.  LAR 
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Pr  =  1.05.  LAR 
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Data  Summary  Tables  -  SAR 


Pr  =  0.44.  SAR 
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Pr  =  1.05.  SAR 
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